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Abstract 

We prove that the multi-time particle distributions for annihilating or coalescing Brow- 
nian motions, under the maximal entrance law on the real line, are extended Pfaffian point 
processes. 

1 Main result. 

Consider a system of annihilating Brownian motions (ABMs) on the real line, where the par- 
ticles move independently except for instantaneous annihilation when they meet. Assume that 
the initial distribution of particles is given by a natural maximal entrance law, which can be 
constructed as the infinite intensity limit of Poisson initial conditions (see [2] or [3] for details). 
The particles, at any fixed time t > 0, form a simple point process on R and it is shown in [3] 
that the (Lebesgue) intensities pt(zi, z%, . . . , z n ) are given by 

p t (zi,z 2 , ...,z n ) = Pf [K t (zi - Zj) : 1 < i,j < n] 

where the Pfaffian is of the 2n x 2n anti-symmetric matrix constructed using the 2x2 matrix 
kernel 

( K}\z) K\\z) \ _ ( -t-Wizt-W) -t-Wizt-W) \ 
Kt[Z) ~ { Kl\z) Kf{z) )~\ t-^F'izt- 1 / 2 ) sgn(z)F(|z|^V2) J 

and F is the Gaussian error function given by 



1 POO 

m= ^l e " 2,ldx 

with sgn(z) = 1 for z > 0, sgn(z) = — 1 for z < and sgn(0) = 0. Briefly, this result was derived 
from a Pfaffian expression for the parity interval probabilities (called spin variables below) for 
ABM's given in H] and [3]. 

This result can also expressed by saying that the positions of ABMs, at any fixed time t > 0, 
form a Pfaffian point process on R with kernel Kf(x — y) (see [5] for an introduction to Pfaffian 
point processes). Using a thinning relation (see [2] or section 2 of [3]) between instantaneously 
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coalescing Brownian motions (CBMs) and ABMs, it is easy to show that under the maximal 
entrance law for CBMs, the particle positions, at a fixed time t > 0, also form a Pfaffian point 
process with the kernel K t (-) replaced by 2K t {-). 

The purpose of this note is to show that the multi-time distributions of particles for both 
ABMs and CBMs can be characterised as extended Pfaffian point processes (see [7J and references 
therein for other examples). 

Notation. We write (G r ) r >o for the semi-group generated by convolution with the Gaussian 
density g r (z) = (2vrr)- 1 /2 e -^/2r > 

Theorem 1 Under the maximal entrance law for annihilating Brownian motions, the particle 
positions at times t > form an extended Pfaffian point process, with multi-time joint intensities 

Pht 2 ...t n ( z ^ z 2,-- -,Zn) = Pf[K(ti,Zi;tj,Zj) : 1 < i,j < n] (1) 
where the space-time kernel K is defined as follows: for t > s and i,j £ {1, 2} 

K ij (t,x;s,y) = G t - a K^(y - x) - 21 {i=1)j=2} g t - s (y - x); 

for t < s and i ^ j € {1, 2}, 

K u (t,x;s,y) = -K u (s,y;t,x), K^(t,x;s,y) = -K^(s,y;t,x); 

and K(t, x; t, y) = Kt{y — x). 

Under the analogous maximal entrance law for CBM's, the particle positions also form an 
extended Pfaffian point process with kernel K(-) replaced by 2K{-). 

Remark. The extra term gt- s (y — x) in the kernel entry K 12 is singular as t J, s, acting like 
a delta function. This reflects the fact that that the particle at (t, y) is likely to have evolved 
from the particle at (s,x). Indeed consider the following heuristic approximation, for small e, 
based on the fact that N s is a simple point measure: 

E[N s {[z,z + e])] « E[(N s {[z,z + e])) 

= lim E[N s ([z,z + e])N t ([z,z + e])] 



rz+e rz+t 

lim / / p st (y,x)dxdy. 

tls Jz Jz 



The left hand side is 0(e) and it is the presence of the delta function that implies the same for 
the right hand side. 

The proof of Theorem [T] is based on the analysis of PDEs solved by certain spin variables 
and particle intensities. We summarize the main steps in the proof in the next section, and refer 
the reader to [3] for certain details. 
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2 Summary of the proof. 



Consider the system of ABM's on R under the maximal entrance law. We write Nt(dx) for 
the empirical measure for the particle positions at time t. Fixing < r\ < < ■ ■ ■ < r m , the 
multi-time intensities ([T]), for t% S {ri, . . . ,r m }, act as intensities for the simple point process 
generated by (N r .(dx) : i = 1, . . . , m) on m disjoint copies of R. See [T] definition 4.2.3 for a 
careful discussion. In particular, for almost all disjoint zi,...,z n we have 



Pt 1 t2...t n {zx,z 2 ,...,z n ) = lime n E 

e4-0 



J\N ti {[z uZi + e]) 



Thus pt 1 t 2 ...t n (zi,z 2 , 



acts a Lebesgue density for the absolutely continuous part of the 



measure -EQliLi ^u(dzi)] (and the measure is non-singular off the set Ui^j{zi = Zj}). To denote 
such an intensity we will use the informal notation 



Ptit 2 ...in( z l> z 2, ■ ■ ■ , z n) — E 



i=l 



Let Nt(0 ■ x) be the number of particles between and x, for any i£R and t > 0. Borrowing 
the terminology from the theory of spin chains [8] we will refer to the following random variables 
as 'spin variables': 

S t (x) = (-l) N t(° :x \ 

In [3] the multi-spin correlation function E[St(xi) . . . Stix^m)} was shown to be given by a Pfaf- 
fian. Taking derivatives in X\, x%, . . . , X2 m -i and then letting a; 2 -i- x Xi ■ ■ ■ > x 2m I ^2m-i leads to 
the intensity (— 2) m p t {xi, x%, . . . ,X2 m -i)- We will follow a fairly similar route, but will include 
an induction over the number of space-time points. 

The following stronger statement, giving a Pfafnan expression for a mixed spin correlation 
and particle intensity, is more suited to an inductive proof. 

Theorem 2 Fix y = (y%, . . . ,y2m) satisfying yi < ])2 < ■ ■ ■ < U2m and times < t\ < t 2 < 
• • • < t n < t, for m,n > 1. Then for ABMs under the maximal entrance law, the intensities 



E 



2m 



3=1 



i=l 



exist, in the sense described in ^jj, and have versions given by the (2n + 2m) by (2n + 2m) 
Pfaffian 

®(ti,zi)(t,y) = {-2) m Pf K(si,Xi;sj,Xj) : (si,Xi),(sj,Xj) G A 
where A is the set of 2m + n space-time points 



{(t,yi), • • • , (t,yzm)} u {(h,z{), (t n ,z n )} 
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and the first 2m rows and columns correspond to the space-time points (t,y±), . . . , (t,y2 m ). The 
kernel K is defined as follows: 

K(ti,Zi;tj,Zj) = K(ti,Zi',tj,Zj) (2x2 entry), 

Kit, yi-,tj,Zj) = (Gt-tjKgizj - yi), G t - tj Kt?(zj - Vi)) (1x2 entry), 

k(t j: z f ,t, Vi ) = {-G t - tj K^{zj - Vi ), -Gt^Kjfizj - yi )) T (2 x 1 entry), 
k(t, yi ;t, yj ) = Kf 2 (yj-yi), for i < j (1 x 1 entry). 

The same result holds if one of m or n is zero, if we take an empty product to have value 1. 

Note that Theorem Q] is the special case of Theorem [2] when m = 0. Furthermore the special 
case when n = was the key to the results in [3]. Note also that the ordering of the entries 
corresponding to the space-time points (tj, z£) in the Pfaffian is not important. Indeed switching 
(U,Zi) for (tj,Zj) will switch two rows and columns at once, leaving the Pfaffian unchanged. 

The main idea of the proof is to examine the equation solved by &r t . :Zj )(t,y) as a function 
of t > t n and y G V2m, where V-jm is the open cell {y : y% < y2 < ■ ■ ■ < y2m}- 

Step 1. Regularity of <£. The regularity of the kernel k implies that (t, y) — > $( t . z .)(t,y) 
defines a bounded function lying in C 1,2 {(t n , oo) x V2 m )nC((t n , oo) x V2m)- The initial condition 
at t = t n may have a jump discontinuity when yj = z% for some i,j. However, the function 
(i,y') -4 $( t . )Z .-)(t,y / ) is continuous as tlt n ,y'^y provided that {yj} je [i t2 m] n {^}i e [i, n ] = 0- 

Step 2. PDE for §u. >Zi \{t,y) satisfies the heat equation 

d 1 

oi®<fii,zi)(t>y) = 2 A$ (*i>*i)(*>y) 

for t > t n and y G V2m- To see this consider expanding $r tiiZj \(t, y) as a finite sum arising from 
the terms of the Pfaffian. We claim each of these terms seperately solves the heat equation. 
Indeed, in each summand, all occurrences of the variables t, (yj) occur inside products of terms 
of the form 

Kf{ Vi - y t ), G t ^K^( Zj - y t ), G^K^zj - y t ). 

Note that each of these terms solves the heat equation. Since each coordinate yi, for i G [1, 2m], 
appears exactly once in each product, the product itself solves the heat equation. 

Step 3. BC for Consider the cell faces F i: 2 m = {y : Hi < ■ ■ ■ < Vi = Ui+i < • • • < y2m} for 
i £ [1, 2m— 1]. On i*2m,i the Pfaffian reduces to a Pfaffian with 2n + 2(m— 1) rows and columns, 
namely where the row and column indexed by (t,yi) and (t, y%+i) are removed. This can be 
seen since these two rows and columns become identical except for the the entries 0,+l,— 1 
where they cross. Then one may subtract row and column (t,yi + i) from row and column (t, y^) 
(using Pf(EBE T ) = Pf(B) det(-E') for the corresponding elementary matrix) and expanding the 
Pfaffian along row (t, yi) leads to the Pfaffian of smaller size. Thus 

*(*,,*,) (f,y) = for y G ^,2™, t > tn 

where y v+1 = (y 1 , . . . , yi~i,y i+ i, y 2 m) (where y hl+l = if m = 2). 
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Step 4. $ as an intensity. Fix smooth compactly supported fa : R — » R, for i E so 
that the distributions (5^ x fa are disjointly supported for i E [l,n]. Consider the integral 



u(t,y) = (f[fa(zi)dz^J $ {t . iZi) (t,y). 



The regularity of 3>(t j|Zi ) above implies that n is a bounded C > ((t n , oo) x V2 m )nC([t n , oo) x T^m) 
solution to the heat equation with initial condition 



(jlMz^dzi^ $(ti,zi){tn,y) for y £ V 2m , 



u(t n ,y) = l_ ( I I <:>,( :,)</.:, ) <!>,/,.,,,(/„. y) Cor y E r 2 ,„. (3) 
and boundary conditions 

«(t,y) = J Rn {f{fa{zi)dz}j $ (t . A) (i,y^ +1 ) for y £ F i)2m , i > t n . (4) 



Step 5. ABM mixed spin correlation and particle intensities. Consider 



v(t,y) = E 



2m 



I[N ti (fa) n%; 

i=i 



1=1 



The simple moment bounds E[\N t ([a, b])\ k ] < C(t, k)\b — a\ k from [3] imply v is a bounded 
function in C([t n , oo) x V2m)- We claim also that v £ C 1,2 ((t n , oo) x V 2 m) and satisfies the heat 
equation on (t n , oo) x Vim- One way to see this is to apply the time-duality from section 2.2 of 
[3] to rewrite 



2 m 



JJ StivAWiN, :s < t n ) 

3 = 1 



E, 



Wm) 



(_l)M^_ tri ,x^_ t ju...u[x t 2 5 t ; 1 ,x t ^ ri 



]) 



where the right hand side is the expectation over an annihilating system of Brownian motions 
started from (yi, . . . , y 2m ) which yields a set of 2K remaining particles at time t — t n , positioned 
at Xf-_ t < . . . < Xt^ tn (here K is random and possibly zero). This satisfies the heat equation 
in (t, y) and can be used to show the same for v. 
Note that v has initial condition 



v(t n ,y) = E 
and boundary conditions 
v(t,y)=E 



2 m 



UN ti (fa) Hs tn ( yj ] 



i=l 



for y £ 



2m i 



(5) 



2 m 



n^(^) n 

1=1 j=l,j^i,i+l 



for y £ Fj >2 m, t > t n . 



(6) 
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Step 6. Induction. We aim to argue inductively that the initial and boundary conditions for 
u and v agree (it is enough to consider boundary conditions only on the faces). This implies 
that u = v and confirms that <&n. tZi \ is the desired intensity, completing the proof. In order to 
use induction we relabel the points as follows. Choose < si < . . . < Sfc and smooth 

compactly supported (fey ■ i G [1, k],i' G [l,nj]), where > 1, so that the distributions 
(6^ x fe : i G [l,n]) are precisely (5 Si X fey '■ i G [1 , fe] , i r G [l,nj]). We argue inductively first in 
k > 0, with a second inner induction on the number 2m of spin space points. Note that when 
k = 0, and the expectation has only spins, the theorem corresponds exactly to the Pfaffian found 
for E[S t (yi) . . . S t (y 2m )] in [3]. 

The initial condition (J5J) for i? can be rewritten, noting that — t n , cis 



i=li'=l i=l 



lim £ 



fe— 1 n. 



2m 



n n ^(<m n^fe) n / m^^^c^c^) 



i=l i'=l 

fc— 1 nj 



i=l 



2~ nk lim £ 



2m 



n n ^(<m n s m n / <4K^u^nO^ 

n n ^i(^,iO n n sm^w 

3=1 



i=l i'=l 



The distributional derivative used in the second equality, (d/dz)St(z) = —2St(z)N t (dz), holds 
almost surely and can be justified as in section 4.3 of [3]. The final expectation can be evaluated 
using the inductive hypothesis in k in terms of the intensities <£. In brief, one can move the 
derivatives back from the test functions fe t i onto these intensities, then carry out the limits 
Z\ \r Zi, and one reaches the corresponding expression for the initial condition ([3]) for u, confirming 
the initial conditions for u and v coincide. More carefully, note that the final expectation 



E 



fc— 1 rii 



2 m 



n n N Sl (fe,i')U s tn(y3)U s MStM) 



3=1 



1=1 



(7) 



is a continuous function of the variables (zi, ■ ■ ■ , z n ,). The integral over R nft can be broken into 
finitely many disjoint regions, according to the ordering of the points (yi : i G [1, 2m]) U (zj, Zj : 
j G [l,nfc]) on the real line. In each such region the inductive hypothesis gives a Pfaffian 
expression for ([7]). These expressions have continuous bounded derivatives in the variables (zj). 
Hence we may integrate by parts to take the derivatives off the test functions 0, , and onto 
the Pfaffian expressions. By the boundedness of the derivatives and the compact support of 
the test functions we may take the limits Zj X z j inside the integral. This leaves a sum of 
Pfaffian expressions indexed over a smaller number of regions, namely the possible orderings of 
the points (y$ : i & [1, 2m]) U (zj : j G [1, n&]). But in each of these regions the Pfaffian expression 
is identical, since switching two of the points yi and Zj, or Zj and Zj, involves interchanging two 
rows and columns leaving the Pfaffian unchanged. Moreover, it is straightforward to check this 
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single remaining Pfaffian is exactly the corresponding expression for the initial condition ([3]) for 
u. 

The above argument shows that the statement of Theorem [5] follows for k and m = from 
the case k — 1. An induction in m implies immediately that the boundary condition ([6]) for v 
agrees with the boundary condition Q for u and completes the double induction. 

Step 7. CBMs. The result for CBMs can be achieved by repeating the above argument 
replacing the spin functional n|=i St(Vj) by the empty interval indicator 

I(N t ([ yi ,y 2 ]) = N t ([y 3 ,y 4 ]) = . . . = N t ([y 2m „ u y 2m ]) = 0) 

throughout. The induction starts due to the Pfaffian shown for the empty interval probability 

Alternatively the thinning relation that connects one-dimensional distributions for ABMs 
and CBMs (under their maximal entrance laws) can be used to check, when y^ < z% for i € [1, n], 
that the spin correlations E [Y\f =1 St^y^St^Zi)] for ABMs coincide with the empty interval 
probabilities E [117=1 ^-(^u([yi, Zj\) = 0)] for CBMs. Differentiating in y{ and then letting Z{ I y^ 
leads to the relation between the intensities (which can be shown to exist via moment bounds). 

Acknowledgement. We are grateful to Neil O'Connell for many illuminating discussions. 
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